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ABSTRACT 

We present an analysis of Chandra observations of two high-redshift clusters of 
galaxies, C1J1113.1-2615 at z = 0.725 and C1J0152.7-1357 at z = 0.833. We find 
C1J1113. 1—2615 to be morphologically relaxed with a temperature of kT = 4.3^g4 keV 
and a mass (within the virial radius) of 4:.3to^ ^ 10^"^ Mq. C1J0152.7-1357, by con- 
trast, is resolved into a northern and southern subcluster, each massive and X-ray 
luminous, in the process of merging. The temperatures of the subclusters are found to 
be 5.5lo!8 5.2lQ;g keV respectively, and we estimate their respective masses to 

be 6.ll};5 X lO^^M© and 5.2^]^;4 x 10^^ Mq within the virial radii. A dynamical analysis 
of the system shows that the subclusters are likely to be gravitationally bound. If the 
subclusters merge they will form a system with a mass similar to that of the Coma clus- 
ter. Two-dimensional modelling of the X-ray surface brightness reveals excess emission 
between the subclusters; suggestive, but not conclusive evidence of a shock front. 

We make a first attempt at measuring the cluster M-T relation at and find 

no evolution in its normalisation, supporting the previous assumption of an unevolving 
M-T relation when constraining cosmological parameters from cluster evolution studies. 
A comparison of the cluster properties with those of nearby systems also finds little or 
no evolution in the L-T relation, the gas fraction-T relation, the /3-T relation or the 
metallicity. These results suggest that, in at least some massive clusters, the hot gas 
was in place, and containing its metals, at z ^ 0.8, and thus that they were assembled 
at redshifts significantly higher than z = 0.8, as predicted in low CIm models. 

We also highlight the need to correct for the degradation of the Chandra ACIS 
low energy quantum efficiency in high-redshift cluster studies when the low energy 
absorption is often assumed to be the Galactic value, rather than measured. 

Subject headings: cosmology: observations - galaxies: clusters: general - galax- 
ies: high-redshift galaxies: clusters: individual: (C1J1113.1-2615, C1J0152.7-1357) 
- intergalactic medium - X-rays: galaxies 
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1. Introduction 

Detailed X-ray studies of distant galaxy clusters have only recently become possible 
thanks to powerful new instruments aboard the Chandra and XMM-Newton satellites. These 
clusters represent the largest bound systems in the universe, and the study of their properties, 
both individually and as a population, allows different cosmological and structure formation 
models to be tested and constrained. 

Many cosmological models describe the formation of structure in the Universe in terms 
of a hierarchical model, in which initial density perturbations (the distribution of which is 
usually assumed to be Gaussian) in the early universe grow under gravity, eventually forming 
the massive structures seen today. In this scenario, less massive systems form first. Then, 
through a combination of merging and the ongoing collapse of the shghtly less dense sur- 
rounding regions, larger systems of galaxies, groups, and clusters of galaxies are formed. The 
immediately observable, luminous components of these systems act as tracers of the under- 
lying distribution of dark matter, which comprises the majority of the mass in the universe. 
While the basic elements of this scenario are widely accepted, many important details, specif- 
ically the epoch and mode of cluster assembly, remain largely unknown. Massive clusters are 
especially useful for quantitative studies because their abundance at high redshifts is very 
sensitive to the matter density of the Universe, their evolution is predicted to be strongest, 
and - being the most luminous systems - they are the ones that are observationally most 
accessible. 

The X-ray properties of galaxy clusters have been exploited as cosmological probes for 
many years, using, for example, their temperatures {e.g. Henry & Arnaud 1991; Henry 2000), 
luminosities {e.g. Jones et al. 1998; Borgani et al. 2001), and more recently, gas mass fractions 
(Allen et al. 2002; Ettori et al. 2003). Since the launch of Chandra and XMM-Newton, these 
properties of clusters can be much more accurately measured out to higher redshifts, enabling 
us to constrain cosmological parameters to much higher precision. Chandra observations in 
particular have also found new and sometimes unexpected features in clusters, such as cold 
fronts {e.g. Markevitch et al. 2000), radio cavities {e.g. McNamara et al. 2001), and shock 
fronts {e.g. Markevitch et al. 2002). The detailed study of individual clusters is an important 
step to compiling samples of high-redshift clusters with well constrained properties. 

Both clusters discussed here, CUIUS. 1—2615 and C1J0152.7— 1357, were discovered in 
the Wide Angle ROSAT Pointed Survey (WARPS: Scharf et al. (1997)); C1J0152. 7-1357 
was also discovered independently in the RDCS (Rosati et al. 1998) and SHARC (Romer 
et al. 2000) surveys. However, due to their high redshifts and the limitations of previous 
instruments, relatively little is known about their properties. Optical spectroscopy was 
performed on galaxies in the region of the ROSAT X-ray sources to measure their redshift 
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and to confirm their cluster status (Ebeling et al. 2000; Perlman et al. 2002). The ROSAT 
observations of C1J0152.7— 1357 have been analysed in some detail (Ebeling et al. 2000), and 
the cluster has also been the subject of both BeppoSAX observations (Delia Ceca et al. 2000) 
and Sunyaev-Zel'dovich effect imaging (Joy et al. 2001). The recent Chandra observations 
of these clusters, presented here, both complement and extend this previous work, and the 
results of the different observations are compared later. These observations, along with 
Chandra and XMM-Newton observations of other high-redshift clusters, will form the basis 
of a study of the evolution of the X-ray temperature function, thereby constraining the value 
of Qm- 

Section 2 describes the general data preparation and analysis methods applied to the 
Chandra data for both objects. Sections 3 and 4 detail the results obtained for C1J1113. 1—2615 
and C1J0152.7— 1357 respectively, before we summarise and discuss these results in section 
5. 

We perform our analysis assuming two different cosmological models: a ACDM model 
with Qm — 0.3, = 0.7 and Hq — 70 km s"^ Mpc"^ , and an Einstein-de Sitter model with 
Hq — 50 km s~^ Mpc~^ and Qm — 1 (^A = 0). The Einstein-de Sitter model was chosen 
to allow comparisons with earlier work. Unless otherwise stated, results are quoted for the 
ACDM cosmology, and all errors are quoted at the 68% confidence level. 

2. Data Preparation and Analysis 

Our two targets were observed with the ACIS-I array of the Chandra X-ray Observatory, 
and the same standard data preparation steps were followed for both clusters as described 
below. C1J0152.7-1357 was observed on 2000 August 8 (ObsID 913) and C1J1113.1-2615 
was observed on 2000 August 13 (ObsID 915). 

Both data sets were reprocessed by the standard pipeline software, removing bad events, 
and excluding ASCA grades 1,5, and 7 (corresponding to diagonal split, "L" -shaped split, 
and 3- pixel horizontal split events). The data were then corrected to the appropriate gain 
map, and for any known aspect offsets. A background lightcurve with time bins of 400 s 
was produced and analysed. The lightcurve of C1J0152.7— 1357 showed only brief peri- 
ods of enhanced background, while there was a significant flaring event in the observation 
of C1J1113. 1—2615. All time intervals during which the background was > 3a from the 
mean level were removed, leaving a useful time of 34 ks for C1J0152.7— 1357 and 88 ks for 
C1J1113. 1—2615. Point sources were identified with the CIAO wavelet detection algorithm, 
wavdetect. 
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2.1. Spectral Analysis 

Spectra were extracted from circular regions centred on the clusters' X-ray centroids, 
within a radius chosen from a preliminary radial profile of the emission to yield good signal- 
to- noise ratios (snr). These radii were typically ~ 50" and are hereafter referred to as the 
spectral radii (r^.) of the clusters. The Redistribution Matrix File (RMF), which relates the 
incident photon energy to the output channel energy of the instrument, and the Ancillary 
Response File (ARF), which describes the dependence of effective area upon energy, both 
vary spatially across the detector. Since our extended sources cover regions described by 
several different RMFs and ARFs, a mean of each of these files, weighted by the source 
counts, was produced and used in the spectral analysis. 

Background spectra were extracted from annuli surrounding the clusters, on the same 
CCD, excluding any point sources and pixels near the chip edges. An alternative method of 
obtaining background spectra from blank-sky datasets (Markevitch 2001) was also examined. 
This method has the advantage that the background spectra are extracted from the same 
region of the CCD as the source, thus eliminating potential systematic errors caused by the 
spatial variations of both the energy response and the effective area across the chip. Changes 
in the quiescent background between the time of the blank-sky observation and that of 
the science observation being analysed can be, approximately, corrected for by normalising 
the blank-sky data appropriately. The method remains, however, vulnerable to temporal 
variations in the spectrum of the particle background, and also cannot easily account for 
the known strong directional variation of the Galactic soft X-ray emission. A comparison 
of both methods showed small (~ la) differences in the values of the best-fit parameters. 
Given the compactness of our targets on an angular scale, we felt that the former method 
of using a background spectrum created from the same observation as the target involved 
approximations that were better understood, and so this method was used in all of the 
analyses detailed below. 

It has recently become apparent that there has been a continuous degradation in the 
ACIS low energy quantum efficiency (QE) since the launch of Chandra. This is believed to be 
due to a build up of hydrocarbons on the optical filter, or the CCDs, which has the effect of 
introducing extra absorption below around 1 keV. In the context of our data, if unmodelled, 
this absorption leads to an overestimate in cluster temperatures. To account for this, an 
XSPEC absorption model ACISABS has been developed^, and the spectra presented here 
were reanalysed with this model. The effect of the QE correction is examined in Appendix 
A. 



^http: / /www.astro. psu.edu/users/chartas/xcontdir/xcont. html 
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The resulting spectra were fit with an absorbed MeKaL plasma emission model (Kaas- 
tra & Mewe 1993; Liedahl et al. 1995) in both Sherpa and XSPEC (Arnaud 1996), with 
completely consistent results. The normalisation of the MeKaL component was used to 
calculate the central density of the intra-cluster gas, pQ, as follows. The MeKaL normalisa- 
tion in XSPEC is defined as 

10-14 f 

" AnDlil + zy J ''-'^HdVcm-^, (1) 

where Da is the angular diameter distance in cm, and rig and uh are the central number 
densities of electrons and hydrogen respectively. Under the assumption that the cluster 
gas has a spherically symmetric density profile described by a /3-profile (Cavalicrc & Fusco- 
Femiano 1976) (the parameters Tc and (3 of which are derived from the data as described in 
the following section), we have 

/ X x2\-3^/2 

p{r) = Po{l+{^] . (2) 



Considering that the spectrum is taken from a circular region of radius r^, viewing the sphere 
in projection, the integral of Eqn. 1 can be expressed as 

Jl.Up^dV = 1-17^^ pirfATir^dr - p{r)HTi (^1 - (^^y^ dr^ . (3) 

Here we assume = l.lTn^. The central gas density po can then be derived from the 
MeKaL normalisation via Eqn. 1. 



2.2. Spatial Analysis 

The spatial analysis of the cluster X-ray emission was performed in the 0.5—5 kcV energy 
range within which the spectra, extracted as detailed in the preceding section, showed to 
feature the highest snr. The spectra were also used to produce spectrally weighted exposure 
maps^ for the imaging analysis. 



^For Chandra data, the standard analysis software produces an exposure map that is an image of the 
effective exposure time in units of s cm^ in the detector plane, taking into account the vignetting of the 

telescope, CCD gaps, the effective area of the telescope, and the dithering of the satellite. By producing a 
spectrally weighted exposure map, in effect combining exposure maps in several energy ranges weighted by 
the number of counts in each range, one also accounts for the energy dependence of the vignetting and of 
the effective area. 
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Two kinds of models of the cluster emission were employed - a one-dimensional (ID) 
model of the radial surface-brightness profile assuming spherical symmetry, and a two- 
dimensional (2D) model of the emission which is fit to an image (see Birkinshaw et al. 
{e.g. 1991) for a discussion of 2D fitting). When producing a radial profile it is desirable to 
have > 20 counts in each radial bin, so that Gaussian errors may be assumed"^. When divid- 
ing a raw-counts image by an exposure map, the latter was normalised to its value at the 
cluster centroid, so that the exposure-corrected image maintains, as much as possible, the 
photon statistics of the recorded counts image. Since our exposure maps vary by typically 
no more than 5% across the region of interest, this approach should not noticeably affect 
the accuracy of our results or their errors. 

A second concern in any imaging analysis is the impact of the point spread function 
(PSF) of the telescope. In our ID analysis we do not account for the PSF at all, because the 
Chandra PSF FWHM is smaller than the typical bin size in the radial profiles, so should have 
little effect. The 2D method deals more elegantly with both the exposure map and the PSF: 
the 2D model is convolved with a PSF model generated for the energy and off-axis position 
of the cluster, and then multiplied by the exposure map, before being compared with the 
original, unfiattened image. This has the advantage of preserving the Poisson statistics of 
the observed image. The disadvantage of the 2D approach is that, since each image pixel 
contains very few (or zero) counts, we must use the C statistic (Cash 1979) to find the best 
fit. Unhke the statistic, the value of the C statistic gives no measure of the absolute 
goodness of fit. The statistic may be used though when fitting a radial profile, provided 
the radial bins contain sufficient counts. Thus the approach followed in our imaging analysis 
was to derive the best fit parameters and their errors with the 2D method, and use ID 
profiles to check consistency and judge the goodness-of-fit. 

In general, we model the surface-brightness distribution of each cluster, or subcluster, 
with a /3-model given by 



where S{r) is the surface brightness at radius r, 5*0 is the central surface brightness, Vc is 
the core radius, and (3 describes the slope of the profile at large radii. In the case of the 2D 
models, the eccentricity (e) and the rotation angle {6) of the semi-major axis of the model 
relate the Cartesian (x,y) coordinates of a point in the image to the radial coordinate. 




3/3+1/2 



S{r) 



(4) 



''For significances below 3a, n > 20 ensures that the Gaussian approximation is accurate to better than 
30 (10) % for the lower and upper limits, respectively, when compared to the exact Poisson errors (Ebeling 
2002). 
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The integrated properties of the cluster gas {i.e. flux, luminosity and mass) were then 
estimated at the virial radius by extrapolating their values, measured within the spectral 
radius, using the measured model profile. The virial radius of a cluster may be estimated 
using self-similar scaling arguments, normalised to numerical simulations. Here, wc adopted 
the formalism of Arnaud et al. (2002a), which may be applied to a flat {Qq + A = 1) or open 
{ilo < 1, A = 0) Universe: 

= 3.804/'A;V2(i + ^)-3/2 
kT N^/^ 



X 



10 keV 



Ko' Mpc. (5) 



Where = (Ac no)/(187r^ fi^), /3t — 1-05 is a normalisation, taken from Evrard et al. 
(1996), and the density contrast Ac is given by 

Ac = 187r2 + 60 (O^ - 1) - 32 {VL^ - if (6) 
for an open Universe, and 

Ac = 187r2 + 82 (O^ - 1) - 39 (O^ - 1)^ (7) 
for a flat Universe. The matter density at redshift z, O^, is given by 

o = ^0 (1 + .o^ 

no{l + zf + {l-no-Iv){l + zf + K- ^> 

The total gravitating mass of the cluster can then be estimated if we assume that the 
cluster is spherically symmetric, in hydrostatic equilibrium, and is well described by an 
isothermal /3-profile. Adopting further a value of 0.59mp for the mean molecular weight of 
the gas, where m-p is the proton mass, we can express the total mass of the cluster within a 
radius r as (Evrard et al. 1996), 

M(r) = 1.13 X 10^^/? J ^ ^''{''"f ^ Mq. (9) 
^ ' keV Mpc 1 + (r/rc)2 " ^ ' 



2.3. Treatment of Errors 

The errors of the measured fluxes (within the extraction regions) were calculated from 
the fractional errors on the normalisation of the spectrum; errors of all derived quantities were 
computed self-consistently with a 'quasi Monte-Carlo method' as follows. For each measured 
parameter T, /3, and flux, an asymmetric Gaussian was computed, centred on the best-fit 
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value and with the positive and negative a widths given by the positive and negative la errors 
found from the best fits. 100,000 vahics were selected at random from these distributions to 
compute the distributions of each dependent quantity, such as Lx or M. The errors quoted 
for each quantity are the range within which ±34% of the values occurred. The core radius 
Tc was fixed during this process. The errors on /3 from the surface brightness fitting were 
computed with Tc free to vary, to account for the correlation between these variables. The 
contributions of the errors on were found to be insignificant compared to those on 6eto, 
and so were ignored. This method highlights the uncertainties involved in extrapolating the 
emission beyond the spectral radius: the extrapolation is strongly dependent on the slope 
of the surface-brightness profile. Lower values of /3 correspond to a shallower profile, with 
a smaller fraction of the emission being directly measured, leading to larger uncertainties 
on the extrapolated values, such as the scaled fluxes. The spread of values of therefore, 
contributes significantly to the errors on the scaled parameters. 



3. Cluster C1J1113.1-2615 

This cluster was identified as a low-surface-brightness, extended X-ray source in the 
WARPS survey (Perlman ct al. 2002). It has a high redshift oi z = 0.725 giving a scale of 
7.2 kpc per arcsecond (ACDM cosmology). In Fig. 1, contours of the X-ray emission detected 
by Chandra arc overlaid on an optical I-band image, obtained at the Keck II 10 m telescope. 
The inner X-ray contours appear smooth and fairly regular, indicating that, at least in the 
central regions, the intra-cluster gas is in a relaxed state. This supports the assumption of 
hydrostatic equilibrium which underpins the derivation of several of the cluster's properties, 
such as the mass and virial radius. The region shown in Fig. 1 contains one faint X-ray 
point source approximately 25" north-east of the cluster. 
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Right ascension 



Fig. 1. — Adaptively smoothed Chandra X-ray contours overlaid on a Keck II (LRIS) I-band 
image of C1J1113. 1—2615. The contours were taken from an exposure corrected image, which 
was adaptively smoothed so that all features are significant at the 99% level. The contours 
start at The contours start at 0.2 counts pixel"^, and are logarithmically spaced by a factor 
of 1.2. 
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Fig. 2 indicates the regions used for the extraction of products. The rectangular regions 
are the source and background regions used in the imaging analysis, while the circle and 
annulus are the source and background regions, respectively, used in the spectral analysis. 
Note that the barred, dashed rectangular region indicates data along the CCD edge that were 
excluded from the spectral background. All point sources were excluded from the following 
spectral and imaging analysis. 
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20 10 11:13:00 

Right ascension 

Fig. 2. — X-ray image of C1J1113. 1—2615 with pixels of 1.968". The regions used for imaging 
(rectangular) and spectral (circle and annulus) analysis are shown. The dashed rectangular 
region was ignored because it was close to a CCD edge. 



3.1. Spectral results 



An initial radial profile suggested a spectral radius of 50", and so a circular region of 
this size, centred on the X-ray centroid {a = 11^13'"5?23, 6 = — 26°15'41 f4) and excluding a 
region of radius 2" around the point source to the north-east, was used to produce a spectrum 
containing around 1000 net counts. The spectrum was grouped to obtain a minimum of 20 
counts per bin, and fit with an absorbed MeKaL model (including the ACISABS absorption 
model discussed in Appendix A), with the absorbing hydrogen column density frozen at the 
Galactic value of 5.47 x 10^*^ cm~^ (Dickey & Lockman 1990). The redshift was frozen at 
0.725, and all fits were performed in the energy range of 0.5 — 8 keV. The spectrum and the 
best-fit model are shown in Fig. 3. 
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The best fit to the spectrum was obtained for a temperature of 4.3lo;4 keV and an 
abundance of 0.62^q;22 solar, with a reduced x^value of 1.00 (62 degrees of freedom). We 
repeated the fit with the abundance frozen at the canonical value of 0.3, and found a tem- 
perature of 4.4tQ5 keV and a reduced x^value of 1.02 (63 degrees of freedom). A redshifted 
iron line is prominent at 3.8 keV. Thawing the redshift parameter of the spectral model led 
to a best fit for a value of z = O.Tibtoxatj which is higher than the optically derived value 
of 2; = 0.725 but not significantly so. Similar redshift differences have been observed in the 
Chandra spectra of other clusters {e.g. Schindler et al. 2001) and are hkely to be due to 
cahbration errors. If so, such errors in the position of the spectral lines will also affect the 
spectral fits of the metaUicities. However, as the differences in the fit parameters caused by 
freezing, or fitting the metal abundance are well within their la errors, we use those values 
derived with the fitted metallicity in our further analysis. 

A hardness ratio map of the cluster was produced by creating images in the energy bands 
0.5 — 1.5 keV (soft band) and 1.5 — 8 keV (hard band), subtracting a constant background 
level in each energy band, and dividing the hard image by the soft image. The images were 
both exposure corrected with the same, broad band exposure map. Strictly, they should be 
corrected by exposure maps generated for the same range of energies as the images, however 
this is a very small effect compared to the large errors due to counting statistics. The low 
number of photons available forced a large pixel size to be used in order to reach > 20 counts 
per pixel, so the resolution of the image was low (110 kpc/pixel). Spectra were simulated, 
using an absorbed MeKaL model at the cluster redshift, for a series of temperatures, and 
the number of photons in the hard and soft bands defined above were divided to give a 
hardness ratio. In this way, a look-up table was constructed, that allowed the conversion 
of the measured hardness ratios into temperatures. The limits of the data allow us only to 
state that the cluster shows no significant departures from a temperature of 5 ± 3 keV on 
scales of ~ 100 kpc. In particular, there is no evidence for a central cooling flow, although 
the constraints are weak. 



3.2. Spatial results 

The cluster X-ray emission was then examined spatially, as described in section 2.2. The 
source model included a constant background flux of 3.95 x lO^^'' counts cm^^ s^^ which was 
estimated from an independent 2D flt to a background region of the same CCD. Four point 
sources in the fltting regions were masked out. After convolution with an appropriate PSF 
model and multiplication with an appropriate exposure map, the model was fit to the data 
in Sherpa using a maximum-likelihood algorithm and the C statistic. The best-fit model 
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parameters were Tc = U.GtH" (105l?6 kpc), P = 0.67tom, and e = 0.2. The background 
level was also allowed to fit, which gave no significant change to the best fitting parameters. 
Errors were computed with all parameter values free to vary, but were not computed for the 
eccentricity because of the computational load involved. 

The image was then divided by a normalised exposure map to produce a exposure 
corrected image (with units of counts), and a radial profile was produced excluding the same 
point sources, estimating the background from the same region as the 2D fitting, and centred 
on the best-fit position of the 2D model. The radial bins were adaptively sized so that their 
minimal width was 2 pixels (ss 1") and the minimal snr value per bin was 3. Emission was 
detected at this level out to a radius of 92.0", which is similar to the size of the region within 
which the 2D fitting was performed. The best-fit parameters of a /^-profile fit to these data 
are = lO.SlJf (76^^ kpc) and /3 = 0.581^;^^, with a reduced value of 1.6 (46 degrees 
of freedom). These fit results do not agree well with those of the 2D fit, however, it is easily 
shown that the disagreement is due to the eccentricity of the cluster emission. An image of 
the best-fitting 2D model, convolved with the PSF and multiplied with the exposure map, 
was used to produce a radial profile in exactly the same manner as above. We obtained best- 
fit values of = 11. 2^];'^" (8llg kpc) and (3 = 0.601qq3, in good agreement with those found 
for the radial profile of the real data. We thus consider the most reliable measurements to be 
those from the 2D fit. The radial profile, best-fitting /3-model, and profile of the best-fitting 
2D model are shown in Fig. 4. 

The reduced statistic of 1.6 indicates that the ID fit is not a very good one, however, 
given that we are fitting a circular model to an elliptical cluster, and that the profile of 
the 2D model appears to agree well with the data, we believe that the 2D model is a good 
description of the data. 
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Fig. 4. — Surface-brightness profile of C1J1113. 1—2615. The sohd hne shows the best-fit ID 
/3-model, and the dashed hne shows a profile of the best-fitting 2D model. 
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To better examine any discrepancies between the data and our /3-model, we subtracted 
the best-fitting 2D model from the data, and smoothed the resulting image with a Gaus- 
sian kernel of fixed size a = 3", thus producing a residual map. The residuals are evenly 
distributed with no significant features, again indicating that the 2D fit provides a good 
description of the data. 

Using Eqn. 5 we estimate the virial radius of the cluster to be 1.3 ± 0.1 Mpc. Our ID 
analysis of the X-ray surface brightness detects cluster emission out to approximately 51% 
of this value, whereas the spectral radius, within which the immediately detected flux is 
measured, corresponds to only ^ 27% of the virial radius. In order to estimate the total 
flux, and hence the total luminosity of the cluster, the observed flux was extrapolated out 
to the virial radius using the best-fit /3-model. To allow comparisons with other work, this 
extrapolation was carried out in both of our chosen cosmologies (the value of is cosmol- 
ogy dependent, as is the conversion between the angular and physical size of the spectral 
radius). We found the scaling factor to be ^ 1.3 for either cosmology. The bolometric, 
unabsorbed, total cluster flux is thus F{r^) = 1.38lo'j2 x 10^^^ erg s~^ cm~^, corresponding 
to a bolometric luminosity of -L(r„) = 3.3lo;3 x 10^^ erg s~^ (ACDM cosmology). 

In order to allow a comparison of the result from the flux measurement with Chandra 
with those of previous observations, we also extrapolated to the virial radius the unabsorbed 
flux in the 0.5 — 2 keV band using our alternative Einstein-de Sitter cosmology. We flnd 
F(r„) = 5.7 ± 0.5 x 10~^^ erg s^^ cm^^, which is marginally lower than the value of 7.4 ± 
1.1 X 10"^"^ erg s~^ cm~^ measured with ROSAT (Perlman et al. 2002). This lower flux 
can be explained by the exclusion of point sources in this analysis that were not resolved 
by ROSAT . Chandra resolves 5 point sources in the region where the ROSAT flux was 
measured. The fluxes of these sources were estimated, assuming the best fltting spectral 
model to the brightest point source within 2 Mpc of the cluster. This was a absorbed power 
law, with a photon index 7 = 2.09 ± 0.05. The combined 0.5 — 2 keV flux in these flve 
point sources was 1.3 x 10"^"^ erg s""*" cm~^, which brings the Chandra flux up to F(r^) = 
7.0 ± 0.5 X 10"^"^ erg s~^ cm~^, in good agreement with the ROSAT flux. 

For further comparison with other work (section 5), the bolometric luminosity was also 
computed for an Einstein-de Sitter cosmology, giving a value of 3.9lo;3 x 10^'^ erg s~^. 

The central hydrogen number density in CUIUS. 1—2615 was found to be uh — 7.3 ± 
0.6 X 10~^ cm~^, giving a gas mass within the spectral radius, and extrapolated out to the 
virial radius, of Mg(r,) = 5.71^;^ x lO^^M© and Mg{ry) = 3.0toi x lO^^M©, respectively. 
The total gravitating mass of the cluster is estimated at M^r^) = 1.11^2 x IO^^Mq within 
the spectral radius of C1J1113. 1—2615; extrapolating to the virial radius we find M(r^,) = 
4.3 ± 0.1 X IO^^Mq. This gives a gas mass fraction (Mg/M) within the spectral radius 
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of 0.05 ± 0.01, and of 0.07 ± 0.01 within the virial radius. A comparison of these cluster 
properties with those measured for a sample of local groups and clusters is performed in 
section 5. 



4. Cluster C1J0152. 7-1357 

Upon its discovery in archival ROSAT PSPC data, C1J0152. 7-1357 was the most X-ray 
luminous distant {z > 0.8) cluster known (Ebeling et al. 2000). Its X-ray morphology as seen 
with ROSAT strongly suggested complex substructure, possibly the result of a merger in 
progress. At the cluster redshift of z = 0.833 1 arcsecond corresponds to 7.6 kpc (ACDM ). 
Fig. 5 shows Chandra X-ray contours overlaid on a Keck II (LRIS) I-band image of the 
cluster. The contours clearly show two peaks of X-ray emission, coincident with two con- 
centrations of galaxies. It appears that the contours are more tightly grouped along the 
south-west edge of the northern subcluster, and more extended to the north-east. This is 
suggestive of the northern subcluster moving toward the southern cluster component. 

It is interesting to note that the inner contours of the southern subcluster are displaced 
to the south of the subcluster galaxies, with the local X-ray centroid ~ 5" south of the 
galaxies. The astrometry of both the optical and X-ray data (using X-ray point sources 
in the field of view) were checked against the Automatic Plate Measuring (APM) source 
catalogue, and found to be accurate to within ~ 0.5"; hence the displacement appears to 
be real. Similar displacements have been observed in other merging galaxy clusters {e.g. 
Markevitch et al. 2002), and can be explained by the coUisionless galaxies (and presumably 
dark matter) moving ahead of the ICM which is slowed by ram pressure. 

The X-ray centroids of the two subclusters (North: a = l'^52"*44?18, S = -13°57'15 ('84; 
South: a = l''52"*39!89, 5 = -13°58'27M8) are 95" apart, equivalent to 722 kpc in our 
ACDM cosmology. The question of whether the two subclusters are gravitationally bound 
is addressed later in this section. 

There is also a region of extended low surface brightness emission to the east of the 
cluster (a = l'^52'"52 ?42, S = — 13°58'5 ''52) associated with an overdensity of galaxies. 
Three galaxy redshifts confirm the existence of a group of galaxies at z = 0.844 ± 0.002 
(Demarco et al. 2002). The limited statistics of both the optical and X-ray data prohibit 
detailed investigation of this system, but it would appear to be an infalling galaxy group. 
Note that this feature is not visible within the field of Fig. 5. 
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Fig. 5. — Adaptively smoothed Chandra X-ray contours overlaid on a Keck II I-band image 
of C1J0152.7— 1357. The contours were produced from an exposure corrected image in the 
energy band 0.5 — 5 keV, which was adaptively smoothed so that all features were significant 
at the 99% level. The contours start at 0.03 counts pixel"^, and are logarithmically spaced, 
by a factor of 1.5. 
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The regions used for the extraction of different products are shown in Fig. 6, overlayed 
on a binned X-ray image. The circular regions were used to extract the source spectra, with 
the elhptical annulus as a background (excluding the dashed, barred rectangular regions). 
The non-barred rectangular regions were used for the imaging analysis. Point sources were 
excluded from all spectral and imaging analysis. 
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Fig. 6. — Chandra X-ray image of C1J0152.7— 1357 with pixels of 2". The regions used for 
imaging and spectral analysis are shown (see text). 
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4.1. Spectral results 

Spectra were extracted within circular regions of radius 49" and 37" of the centroids of the 
northern and southern subcluster respectively, and were fit with an absorbed MeKaL model 
(again including the ACISABS absorption model discussed in Appendix A). The northern 
cluster region contained around 750 net counts, and the southern cluster around 450. The 
absorbing hydrogen column density was frozen at the Galactic value (Dickey & Lockman 
1990) for the position of each of the subcluster centroids (north:1.55 x lO^'^ cm~^, south:1.54 x 
10^° cm~^). The MeKaL redshift was frozen at 0.833 for both chister components, and 
the spectra were fit in the 0.5 — 8 keV energy range, after having being grouped into bins 
containing > 20 counts each. The spectra and best-fit models are shown in Figs. 7 and 8. 
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Fig. 8. — Spectrum and best-fit model of the southern subcluster of C1J0152.7— 1357. 
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The spectral fits yielded a temperature of 5.6l;o;8 keV, and an abundance of 0.15lo!i5 
solar, with a reduced x^value of 0.9 (36 degrees of freedom) for the northern subcluster. The 
southern subcluster was found to be slightly, but not significantly, cooler at a temperature 
of A.S^Qg keV with a metal abundance of 0.701q'44 solar; the reduced x^was 0.7 (35 degrees 
of freedom). We note that the abundances arc poorly constrained and so we prefer to fit for 
the temperature with the abundance frozen at 0.3 solar. This then yields best-fit values of 
5.5^0 7 keV and 5.2t.ol keV for the northern and southern subclusters respectively, without 
significant changes to the reduced values. 

ROSAT observations of this system gave a temperature estimate of 5.9^2'! keV (Ebel- 
ing et al. 2000), while an analysis of a BeppoSAX observation produced a temperature of 
keV and a metaUicity of O.bStoH solar (Delia Ceca et al. 2000). Finally Joy et al. 
(2001) use Sunyaev-Zel'dovich effect imaging of the cluster to estimate its ICM temperature 
at 8.5l^;5 keV. While these quantities were all measured for the cluster as a whole, the X-ray 
temperatures are consistent with that of either subcluster as measured with Chandra . The 
BeppoSAX metallicity is in good agreement with that of the southern cluster and consistent 
with that of the northern subcluster, which is not surprising given the size of the la errors 
of the measurements. The Sunyaev-Zel'dovich temperature estimate is higher than that of 
the Chandra temperature measurement, though at less than 2a. 



4.2. Spatial results 

The spatial distribution of the X-ray emission of C1J0152.7— 1357 was examined as 
described in section 2.2. The system was modelled with two 2D /9-profiles (one for each 
subcluster) and a constant background fiux of 3.37 x 10^^° counts cm~^ s~^ as estimated 
from an independent 2D fit to a background region of the same CCD. Five point sources 
in the fitting regions were masked out, and the model was fit to the data in Sherpa, using 
a maximum-likelihood algorithm and the C statistic. The best-fit parameter values for the 
northern subcluster were = 32.71^^" (248lg kpc), /3 = 0-73tom, and e = 0.08. The fit 
to the southern subcluster yielded Tc = 16.1^2:6" (I22I20 kpc), P = 0.66l:o:o6> and e = 0.00. 
Errors were computed with all parameters free to vary. 

Radial profiles of the X-ray surface brightness were then produced for each subcluster. 
In both cases the profiles were extracted from a semi-circular region spanning 180° away 
from the direction of the merger to avoid contamination by the emission from the other 
subcluster. The radial bins were adaptively sized so that their minimum width was 2 pixels 
(~ 1") and the minimal snr value per bin was 3. In the case of the northern subcluster, 
emission was detected out to 85" at this level; in the case of the southern subcluster, the 
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detection extended to 90". These radii are in good agreement with the size of the 2D fitting 

region. For the northern subcluster a best-fitting radial profile was obtained for Tc = 33^g^" 
{25ltll kpc) and (3 = 0.74toil, with a reduced value of 0.7 (28 de grees of freedom). This 
fit agrees well with the one derived from the best-fitting 2D model; the data and best-fitting 
ID model are shown in Fig. 9. The best-fitting model to the southern profile is parameterised 
by Tc = ll.Qtlf {SStlt, kpc) and f3 = 0.571°;°^, with a reduced value of 0.4 (17 degrees 
of freedom); data and model are shown in Fig. 10. These values are lower than those found 
with the 2D fitting, but not significantly so. In only considering half of the emission from 
each subcluster the errors, already large due to poor photon statistics, are increased further, 
leading to yet greater uncertainties in the ID fit results and low x^values. 2D fitting seems 
to be the best approach to model this complex system as the emission from both subclusters 
can be fit simultaneously, thus making the most of the available photons. This method also 
has the advantage that any excess emission over that expected from two clusters overlapping 
in projection will be apparent in the residuals of the fit. 
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Fig. 9. — X-ray surface brightness profile of the northern subcluster of C1J0152.7— 1357 in 
a semi-circular region spanning 180° away from the direction of the merger. The line shows 
the best-fit ID /3-model. 
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Fig. 10. — X-ray surface brightness profile of the southern subcluster of C1J0152.7— 1357 in 
a semi-circular region spanning 180° away from the direction of the merger. The line shows 
the best-fit ID /3-model. 
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A residual map was produced by subtracting the best-fit model (convolved with the PSF 

and multiplied by the exposure map) from the data, and then smoothing the resulting image 
with a Gaussian kernel of size a = 3". Contours of the positive residuals are shown in Fig. 
11, overlaid on a hardness ratio map of the cluster emission, produced using the regions of 
excess emission as a guide. The hardness ratios of the emission were calculated as described 
in section 3.1. Immediately apparent is a linear feature of excess emission between the two 
subclusters, oriented perpendicular to their merger axis. This narrow feature (width < 5") 
is suggestive of a shock front between the merging subclusters. It is, however, significant at 
only the 4.2(7 level, in a region carefully chosen to maximise its significance (the rectangle 
in Fig. 11; note that this region excludes a point source at the western end of the feature). 
This region also has harder emission than either of the subclusters, corresponding to hotter 
gas (see Table 1), as one would expect in a region of compression, though this result is not 
strongly significant. If confirmed, this feature would constitute the first direct detection of 
a shock in the early stages of a cluster merger. No other systematic residuals are apparent 
in the outer cluster regions, indicating that, overall, our model fits the data well. 

Additional residuals of excess emission are, however, observed in the core of each sub- 
cluster (regions C and F in Fig. 11), suggesting the possible presence of cooling fiows {e.g. 
Fabian 1994), but these residuals are significant at only the 3a level. The core of the northern 
subcluster may have the softest emission in this system (but the errors are large) , consistent 
with a cooling flow, but the opposite is the case in the southern subcluster. Simulations 
predict that cooling flows are disrupted in the late stages of major mergers such as this one 
(Ritchie & Thomas 2002). The possible existence of at least one cool cluster core in this 
developing merger suggests that the central regions of the approaching systems are still fairly 
undisturbed and likely close to hydrostatic equilibrium. 
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Table 1. Values of the hardness ratios, and corresponding temperatures measured for the 

regions shown in Fig. 11 



Region Hardness ratio Temperature ( keV) 



A 


0.82 ±0.08 


o 0+2.6 
O-O-2.0 


B 


0.58 ±0.11 


4.4+}:^ 


C 


1.24 ±0.29 


> 14.1 


D 


0.68 ±0.10 


c, 7+2.0 
' -1.3 


E 


0.78 ±0.14 


7 7+4.8 
'■'-2.6 



-32- 




R^ht ascension 



Fig. 11. — Hardness ratio map of C1J0152.7— 1357, overlayed with contours of excess X-ray 
emission above the best-fitting surface brightness model. Lighter colours represent harder 
emission; the values for each region are given in Table 1. Note that the three point sources 
which appear black were excluded from this analysis 
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Using Eqn. 5 the virial radius of the northern subcluster is estimated to be 1.4±0.1 Mpc. 
The spectral radius used in our analysis thus corresponds to approximately 0.27 r^. For the 
southern cluster we find our spectral radius to reach about 20% of the virial radius of 
= 1.4 ±0.1 Mpc. 

Cluster fluxes for the two cluster components were initially measured within the re- 
spective spectral radii, and then, as before, scaled out to the virial radius by integrating 
the respective best-fit surface brightness model to yield approximate total fiuxes. The unab- 
sorbed total fiuxes were used to estimate the luminosities of the subclusters. For the northern 
subcluster we find a bolometric luminosity of L{ry) = 1.0±0.2 x lO'*''^ erg s^-*^; for the southern 
subcluster the corresponding value is i>(r„) = 5.8^Q;g x 10^'' erg s^^. Combined, these values 
lead to a total bolometric luminosity of C1J0152.7-1357 of L(r^) = 1.6 ± 0.2 x 10^^ erg s'^. 

The subcluster fiuxes measured in the 0.5 — 2 keV band, again extrapolated to r„, 
were combined to estimate the total fiux of C1J0152.7— 1357 in this band for comparison 
with earher measurements, adopting the Einstein-de Sitter cosmology. We find F(r^) = 
l.SStoil X 10-^^ erg s"^ cm-^, slightly lower than the value of 2.9toil x 10"^^ erg s"^ cm^^ 
measured from the ROSAT data (Ebehng et al. 2000). The higher value found in the ROSAT 
data could be explained by a contribution from point sources which were unresolved by 
that instrument, but which are excluded here, and by the fact that the Chandra fiux is 
extrapolated out to the virial radius, rather than to infinity. For completeness' sake we 
also computed the bolometric luminosity of each subcluster in our alternative, Einstein-de 
Sitter cosmology (see Table 3). The total bolometric luminosity of C1J0152.7— 1357 (both 
subclusters combined) is L{ry) = LSllo^gX 10"^^ erg s^^ in this cosmology, in good agreement 
with the value of 2.2^0 J x 10'^^ erg s"^ derived from the BeppoSAX observation (Delia Ceca 
et al. 2000). 

Following the procedure detailed in section 2, the central gas densities, gas masses, 
total gravitating masses, and gas mass fractions of each subcluster were computed, with 
results as hsted in Table 2 and Table 3 for the two different cosmologies. We confirm that 
C1J0152.7— 1357 is indeed a massive cluster; the combined total mass of the two subclusters 
extrapolated to the virial radius is 1.1 ± 0.2 x lO^^M© (ACDM cosmology). 

The Sunyaev-Zel'dovich effect (SZE) work of Joy et al. (2001) quotes an estimate of the 
mass of C1J0152.7-1357 of > 2 x lO^^h^Q^MQ, within a single radius of 65" for the entire 
system. While the different geometry used by Joy et al. (2001) prohibits a direct comparison 
of their estimate with the masses derived here, we note that the masses of the subclusters 
within the spectral radn (~ 50") combine to give 2.4+q;3 x lO^'^h^^MQ in good agreement 
with the SZE estimate. 
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4.3. Cluster dynamics 

We now attempt to address the question of whether C1J0152.7— 1357 is a gravitationally 
bound system. The redshifts of nine member galaxies in each subcluster have been measured 
(Demarco et al. 2002), allowing us to ascribe mean redshifts of z = 0.8352 ± 0.0017 and 
z — 0.8294±0.0013 to the northern and southern subclusters respectively. Via the relativistic 
redshift equation 

^ = (- + ^)^-^ (10) 

the difference between these redshifts can be translated into a hne-of-sight velocity difference 
between the subclusters of Av — 660 ± 258 km s~^. This non- vanishing value of Av can be 
caused by two effects: 

• The subclusters are separated by the difference in the (line-of-sight) distances implied 
by their redshifts, i.e. 58 Mpc, and are being carried apart by the differential of the 
Hubble flow at these distances. 

• The subclusters are at approximately the same distance, and A^; is a peculiar velocity 
reflecting the merger process. 

The former case implies a chance superposition of two separate clusters, which can be re- 
garded as extremely unlikely because the high mass and temperature of the two subclusters 
makes them very rare objects indeed. In the following we assess the probability of the second 
scenario, which implies that the system is gravitationally bound, using the system's observed 
properties and the technique described by Hughes et al. (1996). 

Given a line-of-sight velocity difference Vr and a projected separation Rp between the 
two subclusters of C1J0152.7— 1357, the true pecuhar velocity v and separation R can be 
written as 

Vr — VCOS^y (11) 

Rp = Rsin^R, (12) 

where and are the inclination angles of the velocity and separation vectors, respec- 
tively. Note that different values of and imply an orbital component in addition to 
the radial component of the relative motion of the subclusters. The system is bound when 
the kinetic energy of the subclusters is less than the gravitational potential of the system, a 
condition that can be expressed as 

— sin^RCOs^*^ < 0, (13) 

Rp 
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where M is the total mass of the system. The values of the relevant quantities, for the 
C1J0152.7-1357 system, are Rp = 722 kpc, Vr = 660 km s"^ and M = 1.13 x W^Mq. 
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Fig. 12. — Regions of the (^^,, ^/j) plane corresponding to the bound and unbound solutions 
of Eqn. 13 for the C1J0152.7— 1357 system. The dashed line is equivalent to a purely radial 
motion of the subclusters. 
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Fig. 12 shows the regions of the (^^, "^/r) plane corresponding to the bound and unbound 
solutions of Eqn. 13. The likelihood that the system is bound equals the probability of the 
projection angles having the appropriate values to place the system in the bound region of 
the (\E't,,, \E'i?) plane. This probability is given by the ratio of the solid angle through which 
a given projection angle is observable, to the solid angle of a sphere. Hence the probability 
that has a value in the range ^'i? to + A'^r is given by 



The overall probability that the system is bound can then be obtained as the sum of the 
probabilities of the different combinations of projection angles in the bound region of the 
{"^v, ^i?) plane. Without constraints on the range of permissible angles, we find the proba- 
bility that the system C1J0152.7— 1357 is gravitationally bound, to be 0.78. 



The measured and derived properties of C1J1113. 1—2615 and C1J0152.7— 1357 are sum- 
marised in Table 2, with the properties of C1J0152.7— 1357 given for both the individual 
subclusters and for the cluster as a whole, where appropriate. Cluster properties derived 
using our alternative Einstein-de Sitter cosmology are summarised in Table 3. 




(14) 



5. 



Discussion 



Tabic 2. Summary of the measured and inferred properties of the galaxy clusters C1J1113. 1—2615 and C1J0152.7— 1357, 
based on the Chandra observations discussed here and assuming a ACDM cosmology of ^2^/ = 0.3 {Qa = 0.7) and 
Hq = 70 km s~^ Mpc^^ . Where appropriate, the properties of C1J0152.7— 1357 are listed for both the northern 
(C1J0152.7— 1357N) and southern (C1J0152.7— 1357S) subclusters individually, as well as for the cluster as a whole. The 

temperatures quoted were derived from spectral fits with abimdanccs frozen at 0.3 solar. 



Cluster 


Redshift 


T{keV) 


Lboi(erg s 


rcO 


/3 


e 


rv{Mpc) 


Mg{rv){MQ) 


M{rv){MQ) 


Mg/M{r,) 


Mg/M{rv) 


C1J1113.1-2615 
C1J0152. 7-13578 
C1J0152.7-1357N 


0.725 
0.833 
0.833 


4 O + 0.5 
c 5+0.9 


3.3i;;-;i x lo^* 

5.8-to,l X 10*4 
1.0 ±0.2 X 10"'^ 


32.7111 


fi7+0-03 
"■'"-0.05 

66+° °* 

U.OD_g Qg 

n 70+0.13 
''-'-0.06 


0.2 
0.00 
0.08 


1.3 ±0.1 

1.4 ±0.1 
1.4 ±0.1 


3.0l°;t X 10" 
4.5II ? X 
7.0li;5 X 10" 


4.3I',;;? X 10" 
5.2+J-*^ X 10" 
6.1II J X 10" 


o.o6l°:°? 

0.09l°;°? 


07+0 01 
-0.01 

oq+°o4 

0.12 ±0.04 


C1J0152.7-1357 


0.833 




1.6 ± 0.2 X lO''^ 












1.1 ±0.2 X IQi^ 







''Bolometric luminosity computed from fluxes scaled out to the virial radius. 

I 

00 

I 



Table 3. As for Table 2, but assuming an Einstein-de Sitter cosmology of Hq — 50 km s ^ Mpc ^ and 

= 1 (go = 0.5). 



Cluster F(r„)(0.5 - 2 kcV) erg s-l cm-2 Ltol{erg s'^)'' r^(Mpc) Mg(r^){MQ) M{r^)(M(:-j) Mg/M{rs) Mg/M{r^) 



C1J1113.1- 


-2615 


5.7 ±0.5 X 10- 


-14 


3.9l°;^ X 10^4 


1.1 ±0.1 


^•°-0.3 


X 10i3 


n fj+0.7 


X 10" 


0.06 ± 0.01 


0.08 ± 0.01 


C1J0152.7- 


-13578 




14 


6.7lJ:^ X 10*4 


1.1 ±0.1 


4 0+1° 
^■^-o.s 


X 10^3 


4.3ll:« 


X 10" 


o.o6l°:°? 


0.091°:°^ 


C1J0152.7- 


-1357N 


1.1 ± 0.2 X 10" 


-13 


1.1 ±0.2 X 10"'^ 


1.2 ±0.1 


6-311:? 


X 10" 


5.011:^ 


X 10" 


0.10 ± 0.02 


0.13 ±0.04 


C1J0152.7- 


-1357 


1.8 ±0.2 X 10 


-13 


1.8 ±0.2 X 10''^ 








9-3l?:^ 


X 10" 







''Bolometric luminosity computed from fluxes scaled out to the virial radius. 



-40- 



We find C1J1113. 1—2615 to be a relaxed, hot and massive cluster ai z — 0.725. The 
lack of pronounced substructure suggests that we are observing this system some time after 
its formation, or any recent mergers. On the other hand, the small but non-zero ellipticity 
of the X-ray emission may be interpreted as a remnant of the last merger event. 

C1J0152.7— 1357 is found to be a very massive, hot system ai z — 0.833, which con- 
sists of two massive subclusters that are likely in the process of merging. This scenario is 
supported by our dynamical analysis of the system, which shows that the subclusters are 
likely to be gravitationally bound, the X-ray contours, which show extensions away from the 
hypothesised direction of the merger, and a displacement between the viscous ICM and the 
essentially coUisionless galaxies of the southern subcluster, which suggests a motion of this 
subcluster in a northerly direction toward the northern subcluster. It is difficult to compare 
this observation with recent detailed simulations of cluster mergers as they do not generally 
include galaxies. 

Excess X-ray emission suggestive of a shock front is detected between the subclusters. 
Numerical simulations of equal-mass mergers {e.g. Roettiger et al. 1996; Ritchie & Thomas 
2002) show the formation of such a shock front in the early stages of the merging process 
as the outer regions of the subclusters meet and gas is driven out perpendicular to the 
merger axis. The gas in the shocked region is further compressed and heated as the merger 
continues, becoming detectable by its increased density and temperature. The emission 
detected from the shock region in CIJ0152.7— 1357 is harder than the surrounding emission, 
which is also consistent with the shock front scenario. The overall evidence is suggestive of 
a weak shock front in this merging system, although the limited number of photons prevent 
a robust detection. Future, deeper X-ray observations are needed to clarify the situation. 
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Fig. 13. — Bolometric X-ray luminosities and temperatures of tlic clusters analysed here are 
plotted on the low-redshift L-T relation of Markevitch (1998) in our two cosmologies. Our 
new data points are shown as filled squares (ACDM )and filled triangles (Einstein-de Sitter), 
and represent, from left to right, C1J1113.1-2615, C1J0152.7-1357S, and C1J0152.7-1357N. 



5.1. Evolution of cluster scaling relations 



We now compare the X-ray properties of the three clusters analysed here (we here 
treat the two subclusters of C1J0152.7— 1357 separately) to those expected based on low- 
z scaling relations. We first examine the luminosity-temperature (L-T) relation, using the 
temperatures and luminosities of Markevitch (1998) corrected for cooling flow contamination, 
which should be consistent with the quantities we measure, as the high-redshift systems do 
not appear to host strong cooling flows. The luminosities were derived in our two model 
cosmologies, and are shown in Fig. 13. Previous work, assuming an Einstein-de Sitter 
cosmology has tended to find little, or no evolution in the L-T relation {e.g. Fairley et al. 
2000), and our new data support this; the Einstein-de Sitter points are consistent with 
the local L-T relation. This indicates that contamination by unresolved point sources of 
observations made with earlier satellites (e.g. by Fairley et al. ) was not, on average, a very 
signiflcant effect. Recent work, assuming a ACDM cosmology has found evidence for positive 
evolution of the L-T relation {e.g. Arnaud et al. 2002b; Vikhhnin et al. 2002), i.e. clusters at 
high redshift are more luminous for a given temperature (although Holden et al. (2002) flnd 
no signiflcant evolution in this cosmology). Our data support this general trend of positive 
evolution, albeit with low statistical signiflcance, due to the few points available. 

The loci of our three clusters were also plotted on the low-redshift /3-T relation of 
Sanderson et al. (2002) (Fig. 14), which covers redshifts out to z ~ 0.1. Our three data 
points lie within the scatter around the /5-T relation, suggesting again that the clusters we 
observe at redshifts z ~ 0.8 have a similar distribution of gas to that of systems in the local 
universe. 
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Fig. 14. — (3 values and temperatures of the clusters analysed here plotted on the low-z (3 — T 
relation of Sanderson et al. (2002). Our new data are shown as diamonds, and represent, 
from left to right, C1J1113.1-2615, C1J0152.7-1357S, and C1J0152.7-1357N. The solid line 
shows the best fit to the low-z data, while the dashed line indicates the canonical value of 
/3 = 2/3. 
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The gas mass fractions of all three clusters were measured within their spectral radii, 
and total values were obtained by extrapolating to their virial radii. Note that, since the 
X-ray emission was only detected to 20 — 50% of the virial radii, only around 20% of the 
estimated total gas mass is directly observed. To allow a comparison of our results with 
with those from the low-redshift study of Sanderson et al. (2002) who compute cluster gas 
mass fractions at 0.3r^, we also derived the gas mass fraction of our clusters at this radius 
(slightly larger than our spectral radii) in our ACDM cosmology (this is consistent with the 
values of Hq — 70 km s~^ Mpc"^ assumed by Sanderson and coworkers). The results are 
plotted along with the data of Sanderson et al. (2002) in Fig. 15. The loci of our high- 
redshift clusters do not deviate significantly from the low redshift data, suggesting that, for 
this cosmology, the composition of at least some galaxy clusters does not change significantly 
from z ~ 0.8 to 2; ~ 0.1. 
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0.5 1 5 10 

Temperature (keV) 

Fig. 15. — The loci of the three clusters observed by Chandra discussed here (diamonds) 
plotted on the gas-mass fraction versus temperature plot of Sanderson et al. (2002) for 
low-redshift clusters. For all clusters the gas-mass fraction was evaluated at 0.3r^. The 
new data points represent, from left to right, C1J1113. 1—2615, C1J0152.7— 1357S, and 
C1J0152.7-1357N. 
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5.2. The M-T relation 

Finally, we compare the mass-temperature relation found in the local sample of Sander- 
son et al. (2002) with the properties measured for our high-redshift clusters, and for the 
intermediate redshift sample of Allen et al. (2001). This intermediate redshift sample pro- 
vides a useful comparison, as it has secure mass measurements within i?25oo derived from 
Chandra temperature profiles, and corroborated by weak-lensing measurements. The tem- 
peratures derived by those authors are unaffected by the ACIS low energy QE degradation 
since the low energy absorption was a free parameter in their analysis (Allen et al. 2002). 

Fig. 16 shows the mass (within -R2500) of the hotter (T > 4 keV) systems from the local 
sample plotted against their mass weighted temperatures (measured within 0.3i?2oo)- The 
hue shows the best fitting power law to the local data points, of the form M — ^25007"°, 
with A2500 = 6.2l}:^ X lO^^M^ ^nd a = 1.94 ± 0.12. The intermediate-redshift clusters 
(marked as barred crosses) have their mass (within -R2500) plotted against their mass weighted 
temperatures (also measured within -R250o)- Finally, our new high-redshift data are plotted as 
diamonds, and we plot mass within O.3-R200 against emission-weighted temperatures (under 
our assumption of isothermality, the emission weighted temperature is the same as the mass 
weighted temperature, and constant with radius, though the errors quoted do not include 
any systematic uncertainties in the masses due to this assumption). Note that 0.3it!2oo is 
very close to -R2500 {^-g- Sanderson et al. 2002), so the different scales used are consistent. 

The most important result is that little or no evolution of the M-T relation is observed. 
In order to investigate quantitatively the evolution in the M-T relation we measured the 
mean offsets in mass of the high-redshift points, weighted by their errors, from the local 
M-T relation. The mean offset factors, calculated for the intermediate, and high-redshift 
samples are given in Table 4. The errors quoted are the formal errors in the weighted mean, 
based on the scatter of the data. For both samples the mean offset is consistent with zero, 
a conclusion confirmed by reduced (and corresponding probability) values of 0.40(0.75) 
and 0.46(0.84) for the intermediate and high-redshift samples respectively, when compared 
to the local M-T relation. 

The normalisation of the M-T relation is, however, expected to be lower for objects 
which formed at higher redshifts, because objects which form in a denser universe will be 
hotter for a given mass (T oc M/R). The evolution of the normalisation can be modelled 
by a redshift dependent factor. A cluster of a given mass forming at a redshift z will have a 
temperature given by E{z)M oc T^/^, where E{z) = (1+^)^(1 + zQm + ^a/(1 + - ^a). 
If we assume that the redshift of observation (zobs) is the same as the redshift of formation 
(zf), we would expect the high-redshift clusters to he below the local M-T relation by a 
mean factor 1/E{z) also given in Table 4 (for a ACDM cosmology). In fact, increasing the 
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observed masses by factors of E{z) would place them above the local M-T relation, but 
not very significantly so, given the reduced {and probability) values of 0.62(0.71) and 
3.29(0.02) for the intermediate and high-redshift E{z)M samples, when again compared to 
the local M-T relation (given that the errors on the high-redshift sample are underestimates 
because of the assumption of isothermality) . 

If this lack of evolution in the normalisation of the M-T relation is confirmed by further 
high-redshift cluster observations, then one interpretation is that Zobs 7^-2/, but that clusters 
still reflect the conditions of the epoch at which they formed. In this case the the normali- 
sation of the local M-T relation corresponds to some mean rcdshift of formation of the local 
clusters. The fact that the higher redshift clusters lie on the same, or a similar, M-T relation 
suggests that their redshift of formation is similar to that of the local systems. 

An alternative, and perhaps more realistic, interpretation would be that clusters con- 
tinuously grow in mass, rather than forming in one major merger event identified with a 
formation epoch. In this case the properties of clusters more closely correspond to the epoch 
at which they are observed, and less evolution in the M-T relation is predicted (Voit 2000; 
Mathiesen 2001). A more detailed comparison awaits better data. 
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Tablc 4. Mean weighted fractional offset factors in mass (dM) from the local M-T relation 
of the Sanderson et al. (2002) data above 4 keV, and the predicted offset, 1/E{z), based on 
the assumption of a single redshift of formation for each cluster that is the same as the 
observation redshift. The local offset excludes AWM7, the point far below the relation at 



^ 5.8 keV. 


Redshift Range 


dM 


1/E{z) 


< 0.1 


0.99 ±0.03 




0.1-0.45 


0.95 ±0.05 


0.85 


0.72 - 0.83 


1.03 ±0.08 


0.64 



5 



10 

Mass Weighted Temperature (0.3R2qq) (keV) 



15 



Fig. 16. — The loci of the three clusters analysed here are plotted as diamonds on the low-z 
M — T relation of Sanderson et al. (2002). The solid line is the best fit to the local data at 
T>4 keV and the heavy, barred crosses are from the Chandra data of intermediate-redshift 
clusters of Allen et al. (2001). 
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6. Conclusions 

We have analysed Chandra X-ray observations of two massive high-redshift clusters 
of galaxies: CLJ1113. 1-2615 (z=0.725) and C1J0152. 7-1357 (z=0.833). aJ1113. 1-2615 
appears largely relaxed, whereas C1J0152.7— 1357 is likely to be in the early phase of an 
equal-mass merger of sufficient mass to create the Coma cluster; at 80% confidence the 
system is gravitationally bound. 

Treating C1J0152.7— 1357 as two independent subclusters, we find the hot gas in the 
clusters studied to have very similar global properties to that in massive clusters at low 
redshift. The metallicity, gas fraction and gas density profile slope are all consistent with 
the values found in local clusters of similar temperature, suggesting that the gas was in 
place, and containing its metals, at z — 0.8. These properties point towards an early epoch 
of assembly of at least some massive clusters. 

We measure the M-T relation at z — 0.8 for the first time, and find no significant evolu- 
tion of the normalisation. This behaviour is more consistent with models in which clusters 
are continually evolving, rather than those in which each cluster reflects the conditions at a 
single formation epoch. It also supports the use of the X-ray temperature function at high 
redshifts to constrain cosmological parameters, especially flm. 

Finally, our data support the general consensus of little, or no evolution in the L-T 
relation under the assumption of an Einstcin-de Sitter cosmology, and are consistent with 
recent suggestions of positive evolution in a ACDM cosmology. 
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A. Effect of ACIS QE degradation on spectra 

In this appendix we examine the effect of the ACIS QE degradation due to hydrocarbon 
contamination on the spectra analysed here, and the methods of correcting for this. There 
arc currently three similar methods for correcting for the problem: The ACISABS model for 
XSPEC (http://www.astro.psu.edu/users/chartas/xcontdir/), the same model, implemented 
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for Sherpa (http://cxc.harvard.edu/ciao/threads/sherpa_acisabs/), and corrarf, a Fortran 
code that corrects the ARF using the ACISABS absorption model (http:/ /cxc. harvard.edu/- 
cal/Links/Acis/acis/CaLprods/qeDeg/corrarf.tar.gz). All three methods use the same con- 
tamination rate, which has been derived from calibration observations, and assume the same 
composition of the contaminant. 

We take the spectrum of C1J1113. 1—2615 as a test case, as this is the best quality 
spectrum of the three discussed here. The spectrum was fit with an absorbed MeKaL 
model, uncorrected, and with each of the three correction methods. As the contamination 
increases with time, the correction methods require the number of days between Chandra 
launch and the observation, to normalise the absorption. The observation of C1J1113. 1—2615 
was taken 387 days after launch. The fits were performed in two energy ranges, 0.5 — 8 keV 
and 1 — 8 keV, and the column density was either frozen at the Galactic value, or allowed to 
fit. The best fitting temperatures and column densities are given in Table 5. The spectral 
parameters found for the three correction methods were completely consistent, so we just 
give those for XSPEC ACISABS for clarity. The MeKaL abundance was frozen at 0.3Zq 
for simplicity. 



Table 5. Best fitting spectral parameters with and withont correcting for the Chandra low 
energy QE degradation. Hydrogen columns marked with (f ) were frozen at the galactic 

value during the fit. 



Correction Range ( keV) nH (10^2 cm-^) kT ( keV) Range ( keV) nH (10^2 cm-^) kT ( keV) 



None 


0.5 -J 


i keV 


0.054(f) 


5 2+°-^ 

"^•^-0.6 


1 -J 


i keV 


0.054(f) 


4 7+0.6 

^•'-0.6 


None 


0.5 - f 


i keV 








i keV 




4 7+0-7 
+•'-1.1 


ACISABS 


0.5 - i 


i keV 


0.054(f) 






B keV 


0.054(f) 


4.41°;^ 


ACISABS 


0.5 - J 


S keV 


071+°-°4 


4.31^:? 


1 -J 


S keV 







-53- 



As the QE degradation is most severe below 1 keV we would expect the spectra fitted 
above this energy to give similar results, which is indeed the case. When the fit is extended 
below 1 keV, however, we find the the temperature is overestimated in the uncorrected 
spectra by ~ 15% when the absorbing column is fixed at the Galactic value. On the other 
hand, if the absorbing column is allowed to fit, the temperature found is generally accurate, 
while the absorption is significantly overestimated. It is reassuring, though, that all of the 
corrected spectral fits give consistent temperatures in both energy ranges, and agree well 
with the uncorrected spectra above 1 keV. The best fitting column densities given by the 
corrected spectral fits is also consistent with the Galactic value. 
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